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We investigate thermal photon emission rates in hot hadronic matter from a system consisting 
of TT, p, and ui mesons. The rates are calculated using both relativistic kinetic theory with Born 
diagrams as well as thermal field theory at the two-loop level. This enables us to cross-check our 
calculations and to manage a pole contribution that arises in the Born approximation corresponding 
to the uj —>■ 7r°7 radiative decay. After implementing hadronic form factors to account for Hnite-size 
corrections, we find that the resulting photo-emission rates are comparable to existing results from 
np 7r7 processes in the energy regime of 1-3 GeV. We expect that our new sources will provide a 
non-negligible contribution to the total hadronic rates, thereby enhancing calculated thermal photon 
spectra from heavy-ion collisions, which could improve the description of current direct-photon data 
from experiment. 


I. INTRODUCTION 

One of the goals at the forefront of nuclear physics 
is to understand the phase diagram of strongly inter¬ 
acting matter [1 i- At sufficiently high tempera¬ 
tures, hadronic matter undergoes a transition from color- 
neutral particles to a quark-gluon plasma (QGP) with 
deconfined color charge. As this transition is believed 
to be a crossover, the properties of hadronic matter are 
expected to change continuously as temperatures m- 
proach the pseudo-critical one, Tpc « 160 MeV [1, 1^. 
The phases of QCD matter can be explored in ultra- 
relativistic heavy-ion collisions (URHICs) where a fire¬ 
ball of high-temperature strongly-interacting matter is 
created. This matter rapidly cools and ultimately de¬ 
couples into confined particles which are then detected. 
Hadronic probes of the medium undergo strong rescat¬ 
tering and alterations while traversing the fireball thus 
losing direct information about the properties of the hot 
and dense phases of the evolution. Electromagnetic (EM) 
probes (photons and dileptons), on the other hand, are 
emitted throughout the evolution of the fireball, escaping 
relatively unaltered since their mean free path is much 
longer than the size of the hreball. They thus encode in¬ 
formation on properties of the fireball which are not eas¬ 
ily accessible otherwise, e.g., on its interior temperatures, 
the total space-time volume, and evolution history of its 
collective properties (see, e.g., the reviews in Refs. [1- 
Q for further information on the role of EM probes in 
URHICs). 

Theoretical models of EM probes in URHICs relate 
dilepton and photon spectra to the properties of the emit¬ 
ting medium. The calculated emission spectra of dilep¬ 
tons [^, S agree well with experimental measurements 
thus far [9 [-lllll , while similar calculations of the photon 
spectra [I^ indicate potential discrepancies with recent 
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measurements of direct photons [iM3- This discrep¬ 
ancy concerns both the spectral yields and the photon 
elliptic flow (U 2 ). 

Much theoretical effort has gone into addressing this 
putative “puzzle” [l^, lisl - f^ . The tentative conclu¬ 
sion of most of these approaches is that there are cur¬ 
rently unaccounted-for thermal photon sources from the 
strongly-interacting medium. The difficulty in the prob¬ 
lem and the diversity in possible solutions comes from 
identifying these new sources. For example, additional 
sources from a hot QCP, which contribute early in the 
evolution, are disfavored as they carry a small V 2 (the 
fireball needs a time of the order of the nuclear size, 
TpB Raj to build up the momentum anisotropy) and 
too hard of a spectral slope. This leads one to con¬ 
sider sources which contribute later in the fireball evo¬ 
lution as prime candidates for generating the necessary 
V2 

The present paper aims at identifying such sources by 
investigating novel sources of photon production from 
hot hadronic matter, and in this way contribute to a 
more complete characterization of the electromagnetic 
emissivity of QCD matter in a regime of moderate tem¬ 
peratures. Specifically, we will explore photon emission 
rates from a thermally equilibrated system composed of 
TT, p, and ui mesons, and put their relevance into context 
with existing rate calculations (alternatively our calcu¬ 
lated scattering matrix elements can be used in trans¬ 
port simulations where the momentum distributions are 
not necessarily thermal). The coupling strength of the 
TTpuj vertex, which is pivotal to our analysis, is known 
to be large [ 2 ^ It was instrumental in identifying 
the to t-channel exchange in the np -A Try process as an 
important contributor to the photon emissivity of a hot 
meson gas [3l| (and, in fact, to the w —>■ Tr^y decay M)- 
Thus, a more systematic analysis of pertinent processes, 
including the u! meson as an external particle (incoming 
or outgoing), is warranted. Since the w is not a stable 
particle under strong interactions (it can decay into np or 
Stt), some care has to be taken when evaluating pertinent 
scattering diagrams. To ensure that this is done correctly, 
we calculate the emission rates using both kinetic theory 
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and thermal field theory techniques to identify potential 
issues and cross-check the results. 

The remainder of this paper is organized as follows. In 
Sec.Hl we briefly recapitulate two methods for calculat¬ 
ing thermal photon production in the context of the 'npuj 
system, including the specification of the interaction La- 
grangians. In Sec. nil we calculate photo-emission rates 
using kinetic theory (KT) with Born scattering ampli¬ 
tudes. In Sec. IIVI we use thermal field theory (TFT) to 
compute two-loop diagrams of the photon self-energy to 
both check the KT results and address some intricacies 
in the kinetic approach. In Sec. El we present our results 
and compare them to existing hadronic photo-emission 
rates. We conclude in Sec. ED 


II. THERMAL PHOTON EMISSION RATES 


The thermal photon production rate can be written 
as [13 


dR^ 


-^/(go,T)ImnL(go = WIT) , (1) 


where qo and |g| are the photon’s energy and three- 
momentum magnitude, respectively, and / is the ther¬ 
mal Bose-Einstein distribution function. Employing the 
vector meson dominance (VMD) model [SJ], the 3D- 
transverse part of the EM current-current correlator, 
, is proportional to the in-medium vector meson spec¬ 
tral function (dominated by the p meson) via 


nL( 9 o,<f) = 


(m^°^)^ 


2 ^ ) 
"p 


( 2 ) 


with a dressed p propagator 


-1 


DUqo,q)= qo-W\^ - , (3) 


where is the bare p mass and Ej is the in-medium 
p self-energy. In the VMD model, the thermal photon 
emission calculation thus amounts to calculating the p 
meson self-energy, usually done using TFT. 

Alternatively, calculations of thermal photon emission 
can be done using KT. In Ref. [13, it was shown that 
the imaginary parts of self-energies can be expressed as 
vacuum scattering amplitudes (i.e. using zero-width ex¬ 
ternal particles) folded with appropriate thermal statis¬ 
tical weightings and integrated over the pertinent phase 
space. In this framework, the thermal photon emission 
rate is expressed in terms of al-|-2—>-3-1-7 scattering 
process given by 
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where A/" is the combined spin and isospin degeneracy of 
the incoming particles and |AIp is the initial-state spin- 
and isospin-averaged scattering amplitude of the process 
under consideration. 

The microscopic ingredients to the photo-emission 
rates using the above frameworks are the p self-energy 
(for TFT) and photon-producing scattering amplitudes 
(for KT). Both quantities can be derived on the same 
footing from an underlying interaction Lagrangian. We 
start with a Lagrangian for free n and p fields, 

C^+Cp = ^dp^-d^^-^ml^-^-^pp^-p^'' + \'mlpp-p^ , 

( 5 ) 

with the usual definition of the p held strength tensor as 

Ppv = dpPu - d^pp . (6) 

Interactions between the tt and p mesons are generated 
by minimally coupling derivative terms of the tt and p 
helds to a p gauge held [13, 


dp ^ dp ~\- ippPp ■ T 

dpTT°‘ dpTT°- + 9peabcP|)7r° (7) 

dpPl dppl + gpCabcPlpt , 

where T is the isospin-1 operator with components of 
(T^) be = —Rabc [131 and Pp is the ptttt coupling constant. 
Introducing the above covariant derivatives into the free- 
held Lagrangian results in pirn and ppp interactions of 

TpTZTT — 9pP^ ' {dpTT X tt) , (8) 

Tppp = ~'^9pP^ ' {Pp ^ Pv) ■ (9) 

Interactions with the isosinglet vector lj held are included 
via a Wess-Zumino term [ 23 , [13, 

TwZ = g7vpLjf^^‘'°‘^daUJpdpP^ ■ TT . (10) 

Applying the gauging procedure of Eq. 0 to the Wess- 
Zumino term generates a contact term 

TppTzu) — Pt: pujgp^^ ^de^iOpi^Pp X Pu) • TT . (41) 

Interactions with photons are introduced through 
VMD. We neglect the coupling of the u) held to the elec¬ 
tromagnetic current, as the wy coupling is suppressed by 
a factor of wS-d relative to the py coupling [Ij]. hollow¬ 
ing Ref. , we parametrize the py coupling with Cp, 

Cem = -A^^Cpmlpl , ( 12 ) 

where Ap is the photon held and p° is the neutral p held. 

For a realistic description of interaction processes, we 
also need to account for the hnite size of the mesons, 
which we do in the standard way by introducing pre¬ 
viously determined dipole hadronic form factors at each 











3 


vertex [^1^. For s-channel decay processes, we employ Using the s-channel form factor of Eq. 0, the decay 
at each vertex the form rate becomes 


F{s) = 


2A2 


2 A 2 + [E2{pcm) + E^{pcm)Y 


(13) 


ri-i>2+3 — 


Pcm\MYE'^{pcm) 

Sirnil 


(17) 


where E^Pcm) = y/mf + ^^d pcm{s) is the 

center-of-mass momentum of each hadronic decay par¬ 
ticle, z=2,3, A is an hadronic cut-off parameter (taken to 
be 1 GeV in the present work), and mR is the mass of 
the resonant (or decay) particle. For t-channel scattering 
processes, we employ at each vertex the form 



with t = (pi — p^Y incoming (pi) and outgoing (pY) 
4-momenta (and likewise for w-channel processes). 

Save for s-channel decay processes or those processes 
with individually gauge-invariant vertices, the above im¬ 
plementation of form factors renders it a rather in¬ 
volved task to maintain gauge invariance, especially 
for scattering processes where multiple diagrams con¬ 
tribute [1^. Therefore, we follow the simplified pre¬ 
scription of Ref. by introducing an overall factorized 
form factor squared for each scattering matrix element 
with an appropriately defined average momentum trans¬ 
fer. This is done by identifying the dominant diagram in 
the interaction process, i.e., the diagram with the largest 
contribution to the high-energy photo-emission rate of 
a given process (at low energies the form factor effects 
are small). Typically, this will be a t-channel exchange 
diagram with the lightest exchange particle (e.g., pion ex¬ 
change in TTp 771"), since s-channel processes are sup¬ 
pressed by propagators of the form (s — rnjY)~^- The 
average momentum transfer, i, is evaluated from 


1 


2A2 \ 

2A2 + |t|j 



2A2 \ 

2A^ + \t\) 


(15) 


(and likewise if a u-channel process dominates), and is a 
function of only the photon energy and exchange-particle 
mass in the dominant diagram. This enables us to fac¬ 
torize the form factor from the total amplitude, 

\MY = I Appoint PF(t)4 , (16) 


and thus retain the gauge invariance in the amplitude, 
Appoint j which is evaluated for point-like vertices. 

The coupling constants Qp and g^pu) are evaluated us¬ 
ing data from p ^ tttt and ui 7r°7 decays. In principle, 
the couplings Cp and Pp are related via VMD, such that 
Cp = ejPp. However, for our Born diagrams, where we 
are using zero-width p mesons, we will instead use exper¬ 
imental data from the p —>■ e“'"e“ decay to evaluate Cp. 


where mi is mass of the decaying particle and jAdp is 
the initial-state averaged and final-state summed matrix 
element of the associated decay process. Using Parti¬ 
cle Data Group [1^ values for Fp^^^ = 149.1 MeV, 
= 0.703 MeV, and Fp^g+g- = 7.04 keV, and 
a form factor cutoff of A = 1 GeV, we find coupling 
constants of Pp = 5.98, g^^p^ = 21.6 GeV“^*, and Cp = 
0.0611. 

Having established all necessary interaction vertices, 
including form factor effects and the quantitative eval¬ 
uation of all coupling constants, we now proceed to the 
photo-emission rate calculations. 


III. KINETIC THEORY 

Within the framework of relativistic KT, thermal pho¬ 
ton emission from the 'Kpu system with external lo par¬ 
ticles arises from three 2^-2 scattering processes: Tip 
yw, TTo; —>■ 7 p, and puj jn. Each of the three pro¬ 
cesses is composed of s-, t-, u-channel diagrams and a 
contact (c) term to ensure gauge invariance. The dia¬ 
grams comprising each process in the irpui system are 
shown in Eig. [T] 

The evaluation of photo-emission rates using Eq. o 
requires the coherent sum of the four amplitudes of each 
diagram; 


|Alp — |A4s -l- Ait + A4„ -l- Algp . (18) 

We calculate the amplitudes by applying Eeynman rules 
to the diagrams shown in Eig.[T]and using the Lagrangian 
interactions and form factor procedure detailed in Sec. m 
The explicit expressions for the matrix elements are given 
in Appendix [SI 

Straightforward calculations of rates for the irpuj sys¬ 
tem using KT can be done for the irp yw and puj —yTr 
processes, while the ttoj —>■ yp process reveals a sub¬ 
tlety. In the u-channel diagram, Eig.HKd), the exchanged 
pion can go on-shell, such that u = {pui — P-y)^ = rn^. 
This induces a non-integrable singularity in the phase 
space of the rate calculation using Eq. 0- This pion 
pole configuration corresponds to the uj 7r°y radia¬ 
tive decay, which, in fact, has alread y b een included in 
previous photo-emission calculations jsil . [s^. Therefore 
we must eliminate this contribution from our results to 
avoid double-counting of the radiative uj decay. This is 
facilitated by the structure of the Wess-Zumino term. 


* This slightly differs from the value of 25.8 GeV ^ in Ref. 
since there the VMD value of Cp = e/g = 0.052 is used. 
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(e) 5-channel 


(f) contact term 


(g) t-channel 


(h) u-channel 
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TT 
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(i) s-channel 


(j) contact term (k) t-channel 


(1) n-channel 


FIG. 1: Feynman Born diagrams for photon emission from the npu system, i.e., Figs, (a)-(d) for ttlo 'yp, Figs, (e)-(h) for 
TTp yuj, and Figs, (i)-(l) for pLO —> 'yn. 


Eq. (flUl) . which renders a diagram individually gauge- 
invariant when an outgoing p meson is converted to a 
photon. This feature allows us to separate the u-channel 
from the other three diagrams in the tto; —^ yp pro¬ 
cess while maintaining gauge invariance. Naively, we can 
avoid the w decay by excluding timelike pion configura¬ 
tions with u > 0 from the integration range in Eq. 
but a priori this is not a rigorous justification for this 
choice. To scrutinize this prescription to avoid double¬ 
counting of the ui decay, we turn to TFT where such an 
ambiguity does not occur. 


IV. THERMAL FIELD THEORY 

Thermal held theory provides a rigorous framework 
for the calculation of photo-emission rates. As discussed 
above, within the VMD model the relevant quantity is 
the p meson self-energy; recall Eqs. o and ©• Each 
process and channel considered in the KT calculation has 
a corresponding p self-energy diagram associated with it. 

For our analysis, the relevant self-energy diagrams are 
depicted in Fig. [21 which encompass the u-channel dia¬ 
grams of the pu) 77r and ttw —>■ yp processes. The 
latter is the diagram which involves the divergence and 
double counting of the w —>■ Tr^y decay, while the former 
will be considered to benchmark the equivalence between 
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the two methods. For both cases, we can separate the u- 
channel from the other channels and make a meaningful 
comparison because each diagram is gauge-invariant on 
its own due to the 'npu photo-emission vertex. 




(a) pul —>■ 77r process 



TZ 



(b) TTOJ —>■ p 7 process 

FIG. 2: Cuts of the photon self-energy which generate imag¬ 
inary parts corresponding to the u-channel diagrams of the 
puj —>■ 77r and ttu> —>■ 7 p processes. 


The two diagrams in Fig.[2]have similar structures dif¬ 
fering only by the implementation of the tt self-energy in 
the inner loop. Using standard Feynman rules within the 
Matsubara formalism , the (transverse part of the) p 
self-energy encompassing both diagrams can expressed as 


SJ(9,T) 



xD’'J {uJn,p)D^{qo 


-UJn,q-p) 


where represents the irpuj vertex, p is the three- 

momentum of the w meson, and are the propagators 
of the uj and tt mesons. When evaluating Ej at the 
photon point, the transverse projection operator, , 
can be replaced by the (negative) metric tensor as 

they are equivalent for a gauge invariant self-energy. 

We will focus on calculating the imaginary part of the 
p self-energy, since the real part is small compared to the 
p mass and thereby does not significantly contribute to 
the rates. First, each propagator is expressed in terms of 
a dispersion relation; 


D{po,p) 


1 

TT J-oo Po - W -f le 


( 20 ) 


The (jj is considered a zero-width particle corresponding 
to being an external particle in the KT calculation. Using 


lmD^{po,p) =-nS {pI- p'^ - ml) , (21) 

evaluating the Matsubara sum in Eq. m, and taking 


the p self-energy to the photon point, qo = |(f|, yields 
lmT,p{qo,q,T) = J — 

X |(£’c^go - p- q)"^ IuiDt, {qo - Euu,q-p) 

[\ + r{q^-E^,T) + r{E^,T)\ 

- {Eu,qo + p- g)^Im£)^ {E^ - qo,q-p) 

[r{E^-qo,T)-r{E^,T)]^ , 

where E^ = 

The inner loop constitutes a pion self-energy which fig¬ 
ures in the pion propagators of Eq. (1221) . For interactions 
with thermal mesons, m, this pion self-energy has the 
form [s^ 


V II, f M^m{p,ko,k) 

^TrmV^Oi ^ ) — J 


2En 


X {/’"(£;^,r)-/™(U™ + fco,T)} , (23) 


where M-T^m is the pertinent forward scattering ampli¬ 
tude. Here we consider thermal mesons m = n, p, with 
either tttt scattering through an s-channel p resonance, 
or Tip scattering through an s-channel tt resonance. As 
indicated by the arrows in the diagrams of Fig. [2J taking 
the imaginary part of each p self-energy yields the cor¬ 
responding u-channel Born diagram. The tt self-energy 
Etttt generates the puj —>• yTr process, and the tt self-energy 
Ettp generates the ttlu jp process. 

In the following two subsections we will analyze these 
two processes in more detail. In Sec. IIV A[ we use the 
puj —> ‘jTT process as a benchmark to establish the equiv¬ 
alence of the photo-emission rates calculated using TFT 
to 2-loop order with KT using Born amplitudes. In 
Sec. IIVBI we then use the TFT calculation to identify 
a selection criterion that allows us to eliminate double¬ 
counting with the Lu radiative decay. 


A. pta —>■ 77r a-channel 

While the u-channel diagram of the put ^tt process 
is structurally similar to that of ttlo jp process, the 
key difference is that in the former the exchanged tt can¬ 
not go on shell. This enables straightforward calculations 
for this process using both KT and TFT. The pertinent 
imaginary part of the p self-energy has two major contri¬ 
butions indicated by the two terms in braces in Eq. (l22l) 
and schematically shown in Fig. [31 These two contribu¬ 
tions, commonly referred to as the Landau and unitarity 
cuts (representing p 7 r — >■ w scattering and p — ?► ttlu de¬ 
cay processes for the left-hand-side vertex), can be inter¬ 
preted based on the direction of energy flow of the virtual 
pion. For the right-hand-side vertex, the unitarity cut 
is associated with a pion energy flow oriented into the 
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-^wvw 


a> p y 


FIG. 3: Two cuts of the photon self-energy which give rise to 
imaginary parts corresponding to the u-channel diagram of 
the pu) 77r process. 



qo (GeV) 

FIG. 4: Results from photo-emission calculation of puu —>■ yvr 
via the u-channel diagram at T = 150 MeV. The solid line 
is the result using KT over the full kinematic range of the 
exchanged pion, the dashed line is from TFT via the unitarity 
cut of Fig. [3] and the dot-dashed is from the Landau cut of 
Fig. [3] The sum of the unitarity and Landau cuts is plotted 
but cannot be seen as it coincides with the solid curve. 


irpuj vertex, as shown in Fig.[3ja), and is associated with 
Eui < go- This cut corresponds to first term of Eq. (l2^ . 
The Landau cut is associated with a pion energy flow 
oriented out of the irpu} vertex, shown in Fig. [3{b), and 
corresponds to the second term of Eq. (l22l) . or > qq. 

We have calculated the pertinent imaginary parts of 
the p self-energy given by Fig. [3] and inserted them 
into Eq. to obtain the corresponding photo-emission 
rates; the results for a temperature of T = 150 MeV 
are shown in Fig. |4] and compared to the KT calcula¬ 
tion (both without form factors). We find very good 
agreement between the sum of the Landau and unitar¬ 
ity cuts in TFT and the KT Born amplitude calcula¬ 
tion, thus confirming the equivalence of the two meth¬ 
ods. In addition, we have verified that the Landau (uni¬ 
tarity) cut contribution in the TFT calculation can be 
directly mapped to a calculation in KT where the phase 
space integral is restricted to the energy of the exchanged 
pion flowing out of (in to) the npuj vertex {E^^ > or 
< qo, respectively). The identification of the mapping be¬ 
tween Landau and unitarity cuts and KT for this process 
{puj —> 77r) is facilitated by the fact that the u-channel 
Born diagram does not develop any on-shell singularities 
(contrary to the ttw —)• yp process discussed below). 

Figure [3Kb) shows that the Landau cut gives rise to a 
Born diagram featuring a w —Try radiative decay topol¬ 
ogy. However, the emitted tt is necessarily spacelike as 
it is absorbed by an on-shell p turning into an on-shell 
pion. Likewise, in the Born diagram generated from the 
unitarity cut an on-shell uj absorbs a tt turning into a 
massless photon, requiring a spacelike tt. Similar consid- 


I As in the KT calculations, we use Cp as the EM py coupling in 
Eq. ((TJ instead of the VMD value of ejg. 


erations for the unitarity and Landau cuts will be used 
in the following section to separate out on-shell uj decays 
in the ttlo —?► yp process. 

The TFT calculation affords us with yet another ben¬ 
efit, namely to evaluate the effects of a finite width on 
the p mesons which appear as external particles in the 
Born calculations within kinetic theory. This effect has 
been studied previously within the latter framework and 
found to be negligible [4l|. Within TFT we can render 
this check more rigorous by including a finite width for 
the p meson in the integral for the np loop in Fig. |3Ka). 
We found variations of order 10% or less of the photoe¬ 
mission rates, corroborating the results of Ref. pTj . 


B. TTCJ — >■ yp u-channel 

We now utilize TFT to calculate the photo-emission 
rate given by the p self-energy in Fig. |3Kb) where the 
inner loop is a pion self-energy induced by an interac¬ 
tion with a thermal p meson forming a tt, i.e., from 
Eq. (1331) . Unlike its KT counterpart, this calculation is 
well defined, without any divergences associated with, 
e.g., a pion pole. 

Again, as illustrated in Fig. |5l the p self-energy can be 
separated into unitarity and Landau cuts. However, in 
this process the exchanged tt in the Landau cut can go on- 
shell. When this occurs the Landau cut corresponds to 
the radiative ui decay of w —>■ 7r*^y. This contribution has 
already been included in previous calculations of thermal 
photo-emission rates , and does not constitute a new 
hadronic source from the Trpo; system. 

We now use the distinction between the Landau and 
unitarity cuts provided by TFT as a selection criterion 
for what to include in our rates. We will drop the contri- 
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(a) Unitarity cut 



Tl 

- -► 



P 

—►- 

-^wvw 


a> p y 


FIG. 5: Two cuts of the photon self-energy which give rise to 
imaginary parts corresponding to the u-channel diagram of 
the TTOJ —> 7 p process. 



FIG. 6: Comparison of the contributions to the puj yn pro¬ 
cess from the w-channel diagram (solid line) and the combined 
stc terms (dashed line); no form factors included. 


bution from the Landau cut altogether, corresponding to 
the exchange of timelike pions, to preclude any double- 
counting with the radiative w decay. The remaining con¬ 
tribution from the unitarity cut calculated with TFT is 
found to agree with KT when the integration over phase 
space is restricted to pion energies such that go > E^j . In 
principle, this is a conservative choice since it excludes 
not only all timelike pions but also spacelike ones with 
positive energy. In practice, this difference appears to be 
negligible. 


V. THERMAL PHOTON RATES FROM irpuj 
INTERACTIONS 

Before coming to our final results, let us briefly elabo¬ 
rate on the concrete implementation of form factors for 
the processes at hand. The ttw yp process is dom- 



qa (GeV) 




(GeV) 

FIG. 7: Impact of hadronic form factors on the photo¬ 
emission rates at T = 150 MeV for rrp —>• yuj (upper panel), 
puj —>■ 77r (middle panel), and ttcu —>■ 7 p (lower panel) pro¬ 
cesses. The rates with form factor (solid lines) are compared 
to the ones without form factor (dashed lines). 


inated by the t- and u-channel pion exchanges. Since 
the t- and u-channel form factors have the same struc¬ 
ture, their factorized average form factor is identical and 
is thus applied as overall form factor to all diagrams in 
this process. In the 7rp yto system, the pion t-channel 
exchange is expected to prevail over the u-channel, which 
is suppressed by the p mass in the propagator. There¬ 
fore the pertinent factorized average t-channel form fac¬ 
tor is employed. Before implementation of form factors, 
the poj —> 77r process has two approximately equal con- 


















tributions: the w-channel and the combined contact, s-, 
and t-channel terms (“stc” for brevity), cf. Fig. [51 Note 
that both u and the stc diagrams are individually gauge 
invariant which affords us the possibility to treat these 
two contributions separately. As discussed in Sec. [TTl the 
suppression generated by the averaged factorized form 
factors is driven by the mass of the exchanged particle. 
The u-channel diagram involves an exchanged tt whose 
associated form factor generates a suppression of up to 
a factor 4.5 at Qq = 3.0 GeV. On the other hand, the stc 
term is dominated by t-channel p exchange at high en¬ 
ergies, whose associated average form factor generates a 
suppression of up to a factor 30 at go = 3.0 GeV. Glearly, 
the choice of either from factor would not be an accurate 
procedure. However, due to the separate gauge invari¬ 
ance of the u and stc terms, we can apply an average tt 
exchange form factor (FV) to the u-channel, an average 
p exchange form factor {Fp) to the stc term, and a com¬ 
bination of the two form factors to the interference term 
(which is also gauge invariant), schematically given by 

+ F^F^iMuMtt, + M,tcM:) . 

The net effect of this implementation for the puj —>■ 'jtt 
process is that the total rate is suppressed by a somewhat 
larger magnitude than the other two processes, but still 
less suppressed than if we had used an overall t-channel 
p-exchange form factor. The quantitative effects of the 
form factors are illustrated in Fig. |3 for each of the three 
processes at a temperature of T = 150 MeV. 

One may wonder how variations in the form factor 
cutoff affect the emission rates. The cutoff parameter 
of A^rpuj = 1 GeV (as used throughout this work) has 
been fixed together with the npuj coupling constant in 
Ref. [s^, to reproduce the radiative and hadronic w de¬ 
cays. However, a somewhat smaller cutoff, say, A^rpuj = 
0.8 GeV, still gives agreement with those data once the 
coupling constant is increased accordingly. We have ver¬ 
ified that using the latter set of parameters leads to an 
insignificant change of our photon emission rates over the 
relevant range of photon energies of up to ~5 GeV. 

Our final photo-emission rates for the three processes 
from the npuj system (including form factors) are sum¬ 
marized in Fig.jSjfor three different temperatures of 120, 
150, and 180 MeV. In the phenomenologically important 
regime of photon energies around «1 GeV all three chan¬ 
nels are of comparable magnitude. The puj —>■ 771 channel 
falls off significantly below that, but becomes the dom¬ 
inant photon source at energies above 1.2-1.5 GeV for 
all temperatures. The relative spectral strengths of the 
three channels are rather stable with temperature; only 
the 7 rp —> 7 w rate increases slightly relative to the other 
two at higher energies. We provide parametrizations of 
the final rates in Appendix |B] While the applicability 
of the present methods of thermal photon calculations is 
questionable at temperatures close to the pseudo-critical 
transition temperature, we include parametrizations up 



qo (GeV) 




qo (GeV) 


FIG. 8: Thermal photo-emission rates from the npu system 
for TTp —>■ 'yto (solid lines), ttoi —>■ yp (dashed lines), and pcu —>■ 
yvr (dot-dashed lines) for temperatures T — 120, 150 and 
180 MeV (upper, middle and lower panel, respectively). Form 
factor effects are included. 


to T = 180 MeV for usage in a variety of thermal evolu¬ 
tion models which may have the phase transition imple¬ 
mented at different temperatures. 

At last, we compare the total rate from the newly 
calculated processes with existing literature, specifi¬ 
cally, the u! t-channel exchange in 7 rp ^ 777 as calcu¬ 
lated in Ref. [3l| (Fig. |9l) and tttt Bremsstrahlung from 
Refs. [12,1131 (Fig. (TUI) . The former process involves the 
same vertices as considered in the present work but with 
the ui as an exchange particle rather than an external one; 
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qa (GeV) 



(GeV) 



(GeV) 

FIG. 9: Total rates from the npuj system as calculated in the 
present work (black lines) versus the to t-channel rate (red 
line) for temperatures T = 120, 150 and 180 MeV (upper, 
middle and lower panel, respectively). Dashed lines are with¬ 
out form factor; solid lines are with form factor. 


the pertinent rate was found to be a significant source 
of thermal photons for energies qo > 1.5 GeV relative 
to other known sources in hot and dense hadronic mat¬ 
ter [Sll, H^. Prior to the inclusion of form factors, the 
sum of the newly calculated rates is smaller than the w 
t-channel exchange in the irp jtt process by a factor 
of ca. 4-8 for go > 1 GeV. In the realistic case with form 
factors, however, the two rates are rather close in the 
phenomenologically most relevant range of go < 2 GeV. 
In practice, in URHICs at RHIC and LHC energies, this 


approximately translates into lab-momenta of gt <4 GeV 
due to blue-shift effects of the exploding fireball [l^. In 
Refs. [13,113, tttt Bremsstrahlung was found to be ap¬ 
preciable for photon energies of less than 1 GeV, even ex¬ 
ceeding the contribution from in-medium p mesons with 
baryonic sources [i^, at the lowest photon energies. 
Fig. [TU] shows that the rates from the irpuj system are 
comparable to the Bremsstrahlung rates for photon en¬ 
ergies between 0.5 and 1 GeV, which suggests that the 
contribution of this new photon source is significant rel¬ 
ative to existing thermal photon rate calculations. 

In addition, below the chemical freezeout temperature 
of Tch — 160 in URHICs, effective meson chemical po¬ 
tentials build up, in particular for pions, which further 
augment two of the three newly calculated rates. The 
TTU! -P- 7 p and puj -P- "fir processes will pick up pion fu- 
gacity factors Zt^ = exp(/i^/T) to the 4*^ and 5*^ powers, 
respectively, compared to the 3'^'^ power for irp -p- jtt pro¬ 
cesses. With this enhancement from pion fugacities, the 
thermal photon emission processes calculated could pro¬ 
vide a non-negligible contribution to the direct-photon 
spectra in URHICs [^ . 

In our current study, we have focused on the npuj sys¬ 
tem due to its known large coupling constant and the 
relatively small masses of these hadrons. Decay con¬ 
tributions of higher resonances as intermediate particles 
in hadron-p scattering processes have been considered in 
Ref. [13 and found to give subleading contributions at 
the photon point (cf. Fig. 1 in that reference). The ques¬ 
tion remains whether the scattering of higher-mass states 
can give significant contributions. As an example, let us 
consider the oi as an external particle. This is the next 
higher-mass particle with a known large coupling to np 
and Try, and closely resembles the npuj system in the per- 
tinent photon-emitting processes. Inspection of Fig. I of 
Ref. [ 34 I reveals that the oi contribution to the p self¬ 
energy is down by about an order of magnitude compared 
to the u! at the photon point. This factor can be readily 



qo (GeV) 

FIG. 10: Total rates at T = 150 MeV from the Tvpuj system 
(solid black line) c omp ared to the tttt Bremsstrahlung rate 
(dashed black line) |42.l43 and the rrp yrr uj t-channel rate 
(red line). 
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understood by realizing that the Trpw coupling is a factor 
2 larger than the irpai coupling, and that the thermal 
oi density is suppressed by more than a factor of 3 at 
temperatures of T = 150 MeV. 

VI. CONCLUSIONS 

In this work, we have calculated the photo-emission 
rates from the tree-level scattering processes irp —> yw, 
TTw —>■ 7 p, and poj —>■ yTr using relativistic kinetic the¬ 
ory. Complementary calculations were performed using 
thermal field theory for the w-channel diagrams of two of 
the processes. This allowed us to (a) explicitly establish 
consistency between the two methods (as well as exert 
quality control of the results) and (b) identify a crite¬ 
rion by which to exclude singular contributions from the 
exchange of on-shell pions in ttw —> yp and thus avoid 
double-counting with the w radiative decay. After the 
inclusion of hadronic vertex form factors, which suppress 


the rates at high energies, our total rate resulting from 
all three processes turns out to be comparable to the ui 
t-channel rates in the irp —> yyr process and to tttt 
B remsstrahlung [d^. Our results thus provide an en¬ 
hancement of existing rate calculations in the direction 
of what has been conjectured in Ref. (^ : we anticipate 
that thermal emission from irpuj interactions will give a 
sizable contribution to direct-photon spectra in heavy-ion 
collisions. The precise extent to which the new thermal 
photon sources may help to reduce discrepancies with 
experimental spectra and elliptic flow will require their 
implementation into evolution models for the fireball in 
heavy-ion collisions. Work in this direction is in progress. 
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Appendix A: Born Matrix Elements 


In this appendix, we write out the matrix elements corresponding to each Feynman diagram in the irpuj 
system, explicitly indicating isospin indices and labeling each particle’s four-momentum. They are derived by ap¬ 
plying standard Feynman rules to the diagrams shown in Fig.[T]and using the interaction Lagrangians of Eqs. (l8l)- (IT^ . 

• 7 r“(pi) p^{p 2 ) 7 ( 9 ) process with isospin indices 7 r“p^: 


7 J'“(pi) 7(9) 



FIG. 11: Diagram showing four-momenta labels and isospin indices for the process np —>■ yo;. 


Ms = {2pi -P2)sq^q°‘ e^ja/3 {p2)e*'^{q)e*^(ps) 

s — 

Mt = [2pi - 9)tP^9“ £\p2)e*''{p^) 

i 

- inn - \ 

t — \ J 

[-957(^2 + q)x - 9 j\{p 2 - 29)^ + 95 a( 2 P 2 - 9)7] (P2)£*'^iq)£*^(Ps) 

Me = -9TTpuj9pCpq°'e^5ape\p2)£*"*[q)£*'^{P3)(?°‘^ 


(Al) 

(A2) 

(A3) 

(A4) 


• 7r“(pi) a;(p 2 ) ^ liq) P^iPa) process with isospin indices n°‘p^: 
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7r“(pi) 7(g) 



FIG. 12: Diagram showing four-momenta labels and isospin indices for the process ttw —>■ 7 p. 
^ (g + P.)-(g + P3)n 

5 — \ / 

[5A5(g + 2p3)j - g\j{2q + P3)s + g.ys{q - P3)\] e'^(p2)e*'^(g)e*‘^(P3)e^“*’ 

_ g)Pp“(2pi - g)..^ e^5a/3 {p2)e*'^{q)e*^{ps) 6^“'’ 

C TO.^ 

- g)>2 (pi - P2 + g)^ {P2)e*"'{q)£*^ {P3) 

A4c = g^pujgpCpP2£5^ap£^{jP2)£*'^{q)£*^{p3)£^°‘'" 


(A5) 

(A6) 

(A7) 

(A8) 


• P°‘{pi)^{P 2 ) —>■ 7(g) T^^ips) process with isospin indices p“7r^: 

p“(pi) 7(g) 



FIG. 13: Diagram showing four-momenta labels and isospin indices for the process puj ^ 'pn. 


Ms 


Mt 


Mu 

Me 


g.p^gpCp ^p^^^^ ^ 2P3)^ £\pi)£^ip2)£*^iq) 

s 

- ,)>; (-,•'>■ + (njiiOBjiA) 

t — TOp \ rrip ) 

[-5<57(f’i + 9)^ - g-i\{Pi - 2g)5 + 55 a(2pi - g)^,] tpvap £\pi)£^ip 2 )£*'^(q) 

pi^p^{p^ - q+p3)s e p e'^(pi)£^(p2)e*^(g) e^“'' 

u — 

g-KpuigpCpP2 esjap e'^(pi)e^(p2)e*'^(g) e^“'' 


(A9) 

(AlO) 


(All) 

(A12) 


Appendix B: Parametrizations 

In this appendix, we present parametrizations of the photo-emission rates for each process in the irpu} system, 
along with plots of comparisons of parametrizations to calculated rates. We have verified the accuracy of the 
parametrizations to within 10% for temperature and photon energy ranges of 100 MeV < T < 180 MeV and 
0.2 GeV < go < 5.0 GeV, except for the lowest photon energies of the pui —>■ 77r process, whose overall contribution in 
that photon energy range is negligible. Form factor effects are included in all rate parametrizations. 
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Tip —)■ 70 ; 



?o (GeV) 



FIG. 14: Left panel: Calculated thermal photo-emission rates from the 7 rp ^ 70 ; process (symbols) compared to the 
parametrized rates (lines). Right panel: Ratio of parametrized rates to calculated rates. 


ai(T) 

«2(T) 

«3(T) 

ai{T) 

05 (T) 
aeiT) 
ar{T) 


exp [ai(?o-I-02 + a3(?Q^-I-05(90 + 05)“^] [fm GeV , 

-35.8991 -h 460.425 T - 2592.04 + 5342.32 , 

-41.9725 -h 601.952 T - 3587.8 -h 7604.97 , 

0.740436 - 16.7159 T + 133.526 - 347.589 , 

2.00611 - 3.79343 T + 29.3101 - 72.8725 , 

-8.33046 -h 121.091 T - 801.676 + 1712.16 , 

17.9029 - 388.5 T + 2779.03 - 6448.4 , 

-15.622 -h 340.651 T - 2483.18 + 5870.61 . 


(Bl) 


(B2) 


• TTUJ —)■ 7P 




FIG. 15: Left panel: Calculated thermal photo-emission rates from the tiuj —>■ 7 p process (symbols) compared to the 
parametrized rates (lines). Right panel: Ratio of parametrized rates to calculated rates. 
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90 


'TTLO—^-yp 

IFq 

ai(T) 

«2(T) 

a3(T) 

04 (T) 
a5(T) 
06 (T) 

ar(T) 


exp [oigo + 02 + Oagp* + 05(90 + 06)“’'] [fm GeV , 

-29.4663 + 291.356 T - 1301.27 + 2102.12 , 

-45.081 + 688.929 T - 4150.15 + 8890.76 , 

-0.260076 + 8.92875 T - 60.868 + 136.57 , 

2.2663 - 8.30596 T + 49.3342 - 90.8501 , 

10.2955 - 317.077 T+ 2412.15 - 6020.9 , 

3.12251 - 47.5277r + 222.61 - 241.9 , 

-3.39045 + 56.5927 T - 336.97 + 622.756 . 


(B3) 


(B4) 


• PUI ^ 'fTT 


> 

a 

O 

O ' 

n 

3 

cl' 

■a 



0.2 0.5 


1.0 


1.5 


2.0 


2.5 


3.0 



qo (GeV) 


qo (GeV) 


FIG. 16: Left panel: Calculated thermal photo-emission rates from the pu —>■ 'yn process (symbols) compared to the 
parametrized rates (lines). Right panel: Ratio of parametrized rates to calculated rates. 


dRpu)^'^-K 

ai(T) 
02 (T) 

03 (T) 

04 (T) 


exp [ar9o + o, + GeV^] 

-29.6866 -h 331.769 T - 1618.66 -h 2918.53 , 

-15.3332 -h 90.2225 T - 300.185 -h 428.386 , 

-7.35061 -h 109.288 T - 630.396 -h 1227.69 , 

-10.6044 -h 109.1 T - 500.718 + 872.951 . 


(B5) 


(B6) 
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